function [zt,theta] = cpoints(ntotal,r)

% Generation of the intial points over the curve zshape.
% The points are (angular) uniformly distributed.
% Even generation for circular shape, the angle between -pi and pi.
 
 dth=2.d0*pi/ntotal;

      for ii=1:ntotal
         theta=-pi+(ii-1)*dth;         
         zt(ii)=zshape(r,theta);
      end

      end
